1/f Flux Noise in Josephson Phase Qubits 
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We present a new method to measure 1/f noise in Josephson quantum bits (qubits) that yields 
low-frequency spectra below 1 Hz. Comparison of noise taken at positive and negative bias of a 
phase qubit shows the dominant noise source to be flux noise and not junction critical-current noise, 
with a magnitude similar to that measured previously in other systems. Theoretical calculations 
show that the level of flux noise is not compatible with the standard model of noise from two-level 
state defects in the surface oxides of the films. 

PACS numbers: 



Superconducting integrated circuits are a leading 
candidate for scalable quantum information process- 
ing (QIP) PQ. Quantum bits (qubits) based on Joseph- 
son junctions have already achieved several key mile- 
stones, including single and coupled qubit state tomog- 
raphy [H[3]. Moreover, the dominant mechanism for en- 
ergy relaxation is becoming understood @], and steady 
improvements can be expected in the coming years. How- 
ever, to realize the full potential of Josephson junctions 
for QIP, it will be necessary to extend qubit dephasing 
times. Present dephasing times are in the 100's of ns 
range; the short coherence places a strict limit on the 
number of gate operations which can be implemented, 
and represents a significant obstacle to scaling up. De- 
phasing is produced by low-frequency fluctuations in the 
qubit energy. In the case of the Josepshon flux qubit and 
the flux-biased Josephson phase qubit, these fluctuations 
are believed to arise from a magnetic flux noise applied 
to the qubit loop, with a spectral density that scales in- 
versely with frequency {1/f). Moreover, the magnitude 
of the flux noise inferred from qubit Ramsey fringe exper- 
iments is of the order of several zx<I> / -\/Hz for both three- 
junction flux qubits and phase qubits, despite a difference 
in loop inductance of almost two orders of magnitude [5]. 

Low-frequency noise in superconducting circuits has 
been studied for decades in the context of amplifiers 
based on the Superconducting QUantum Interference De- 
vices (SQUIDs) [6j|7]. More that 20 years ago in a series 
of experiments on SQUIDs cooled to millikelvin tempera- 
tures, researchers found that the devices displayed a flux 
noise with a power spectrum which scaled like l/f a at 
low frequencies, where a lies in the range from 0.6 to 1. 
The magnitude of the noise was seen to be only weakly 
dependent on a wide range of device parameters such as 
SQUID loop inductance, geometry, material, etc., with 
a canonical value at 1 Hz of about 2 /i$o/\/Hz. The 
origin of the excess low temperature flux noise in these 
experiments was never understood, and the issue has lain 
dormant for almost two decades. Now it seems clear that 
the excess low-temperature noise of these SQUIDs is in- 



timately connected to the measured dephasing times of 
superconducting qubits [SJ [51 [S] • 

In this Letter we present the results of a novel measure- 
ment in a Josephson phase qubit that uses the resonant 
response of the qubit to directly measure the spectrum of 
low-frequency noise. This general method can be used for 
any qubit system. By alternating the sense of the qubit 
bias, we show that the noise is predominantly flux-like, 
as opposed to a critical-current noise. This experiment is 
the first to directly connect flux noise in superconducting 
qubits to previous measurements in SQUID devices. Ad- 
ditionally, we present the results of calculations of flux 
noise from paramagnetic defects in the native oxides of 
the superconductors, and show that the measured flux 
noise is not compatible with the standard model of two- 
level state (TLS) defects. 

A photomicrograph of our device is shown in Fig. [T] 




FIG. 1: Photomicrograph of Al-based qubit fabricated on a 
sapphire substrate using a SiN^ dielectric for crossover wiring. 
(A) Josephson junction with area Aj ~ 2 /im 2 and critical 
current Iq — 1.9 fiA. (B) Qubit inductor with inductance 
L = 800 pH. (C) Readout SQUID. (D) Qubit flux bias. (E) 
Qubit shunt capacitor with C = 1 pF. (inset) Qubit potential 
energy U as a function of the superconducting phase S across 
the qubit Josephson junction. 
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FIG. 2: (A) Qubit response curve for the probability Pi of the qubit in the |1) state versus microwave excitation frequency at 
the positive flux bias $+. (B) The time evolution of Pi is measured for at the two frequencies ojJ and uj, which are shown 
in (A). (C) and (D) are the same as (A) and (B), respectively, but are for In (B) and (D), data was taken at approximately 
800 samples per second. The correlation in uy£ and ujZ indicates flux noise. 



a more detailed discussion of its operation is given else- 
where .10. A current bias I = §/L is applied to the 
Josephson junction via a flux $ threading an inductor 
L placed across the junction. The bias current is set 
slightly below the critical current of the junction I so 
that the system can be well modeled by a cubic potential 
(see inset). The two lowest quantum states in this poten- 
tial well are labeled as qubit states |0) and |1), and have 
an energy difference E\q that can be tuned with bias. 
Transitions between |0) and |1) are driven by applying 
microwaves at a frequency loio/2tt — Eio/h ~ 6 GHz. 
The qubit state is measured by applying a fast bias pulse 
to lower the potential barrier, forcing only the |1) state 
to tunnel out of the well [11] . 

The transition frequency wio is given by 



wio ^ uJpo (1 - \I\/Io) 



1/4 



(1) 



where w p o — 2 1 / 4 (27r/ /C<I>o) 1 / 2 . Low frequency fluctua- 
tions in the current bias / and critical current Iq produce 
fluctuations in the transition frequency primarily from 
the second term. The qubit can be operated at both 
positive and negative current bias. A positive fluctua- 
tion in I gives an increase in lu w at both positive and 
negative current bias, a symmetric change. A fluctuation 
in /, however, gives an asymmetric change. Therefore, a 
spectroscopic measurement of the transition frequency 
at positive and negative bias currents provides a clear 
differentiation between these two different noise sources. 
Moreover, the steep response of the resonance allows for a 
reasonably sensitive measurement of the fluctuation mag- 
nitude. 

The experiment is performed by choosing positive and 



negative current biases close to the critical current, cor- 
responding to flux biases and that have approx- 
imately equal transition frequencies. A spectroscopic 
measurement is then performed by applying a long 2 ^s 
microwave pulse and measuring the probability P\ of the 
occupation of the |1) state. The amplitude of the mi- 
crowave excitation is chosen so that Pi < 0.4 at peak 
response to prevent significant power broadening of the 
qubit response. Qubit response curves for positive and 
negative biases are shown in Figs. |2}\ and[2p. 

As shown by Eq. ([TJ, fluctuations in bias and critical 
current will cause these resonance curves to shift. The 
probability P\ is most sensitive to qubit bias at the half 
maximum points of the resonance curves, labeled as fre- 
quencies u>^, uj^ and wj, uj^ in Figs. [2j\ and re- 
spectively. A plot of Pi versus time is shown for these 
four frequencies in Figs. [2j3 and[2|3. An anti-correlated 
change in Pi within the data pairs (u>2, w^) and (a/jr, 
a/ft) is expected, and represents a systematic check of 
the measurement method. The small deviations from 
anti-correlation are due to other influences, such as fluc- 
tuations of resonant TLS defects [IT] which affect the 
measurement probability pairs in a correlated manner. 

The data at and u>~^ give a symmetric correlation of 
Pi with time, as can be seen from the traces in Figs. [2j3 
and [2p. This shows that the dominant low- frequency 
noise for the qubit is a flux noise. The relation between 
Pi and flux is calibrated by measuring Pi while sweep- 
ing the qubit flux bias for each of the four frequencies. 
The flux noise data measured at u>l and lur are aver- 
aged for positive and negative bias, then Fourier trans- 
formed, cross-correlated, and averaged over frequency to 
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obtain the cross-correlated flux noise spectrum plotted 
in Fig. |3}\. We note that white noise from the measure- 
ment process is automatically subtracted in this cross- 
correlation analysis. The noise has a 1//™ spectrum with 
a = 0.95 and extrapolates to a flux noise at 1 Hz of 
4 /i$o/\/Hz- This magnitude is comparable with previ- 
ous measurements oil/ f flux noise in superconducting 
devices [6j|7]. 

In Fig. [3]3 we plot the correlation amplitude and 
phase angle of the cross-correlated flux noise spectrum. 
The correlation angle of zero indicates asymmetric (flux- 
like) noise. At the lowest frequencies, the contribu- 
tion from measurement noise is small. The satura- 
tion of the correlation amplitude at a value slightly less 
than unity may indicate a small contribution from junc- 
tion critical-current noise. Taking the contribution from 
critical-current noise to be ~ 5%, we find S/ (l Hz) = 
0.05 Hz)/L 2 = 1.4xl(r 12 $/Hz for our 2 ^m 2 area 
junction. This value is compatible with previous exper- 
iments [7J, giving a critical-current noise at 20 mK that 
is a factor of 36 lower than predicted for 4.2 K [T2"] . 

We note that flux noise produces dephasing of the 
qubit state, as can be measured directly in a Ramsey 
fringe experiment. The magnitude of our noise is within 
a factor of two of that required to explain our qubit de- 
phasing times of around 200 ns [8] . 

In what follows, we examine the possibility that the 
flux noise is due to magnetic TLS defects in the native 
oxides of the superconducting films, as was recently pro- 
posed in Ref. [9]. The standard TLS model [13] de- 
scribes an ensemble of defects, each with two micro- 
scopic configurational states \L) and \R) that have a 
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FIG. 3: (A) Cross-correlated noise power spectrum of data 
taken from Figs. [2|B) and (D). The line for an ideal 1// 
spectrum is shown for reference. (B) Correlation amplitude 
(diamonds) and phase angle (squares) of cross-correlated data 
in (A). 



two-state Hamiltonian with diagonal matrix elements 
±A/2 and off-diagonal elements Aq/2 due to tunnel- 
ing. The eigenstates are given by \g) = sin(f?/2)|£) + 
cos(0/2)|i?) and |e) = cos(0/2)|L) - sin(0/2)|i?), where 
9 = arctan(Ao/A). The difference in energy of the 
two states is E = \/A 2 + A§. The defects are as- 
sumed to have a constant distribution of energies A, 
but a log-uniform distribution in Ao because tunneling 
is exponentially dependent on parameters. Upon chang- 
ing variables to (E, sin0), the joint distribution is given 
by d 2 N — p dE dsin0/sin0cos0, where p is a ma- 
terials constant describing the defect density of states. 
Dipole radiation of the TLS via phonons gives a relax- 
ation rate determined by the matrix element sin0, yield- 
ing Ti = r™ ax sin 2 8. The resulting log-uniform distribu- 
tion of Ti produces a 1// noise spectrum. 

To estimate the magnitude of the flux noise from mag- 
netic TLS defects, we consider a TLS magnetic moment 
equal to the Bohr magneton /i b , and further assume that 
fluctuation of the TLS will completely randomize this 
magnetic moment (we discuss the validity of this assump- 
tion below) . Following the analysis of Ref. [T3] , one can 
show that the low-frequency spectral density of the TLS 
magnetic moment per unit volume is given by 



S m (uj/2n) ~ 4kTfj? B p 

_ kTp 2 B p 
~ w/2tt 



<2Ti 2Tx 
2Y[T 2 +oj 2 



(2) 
(3) 



In order to connect the above expression to the mea- 
sured flux noise, we need to know how each TLS couples 
magnetically to the SQUID. Analytical expressions for 
flux noise may be calculated using reciprocity: the mag- 
netic flux from a spin of moment m is given by (B ■ m) /I, 
where B is the magnetic field at the spin produced by a 
test current / in the SQUID loop. We consider two ideal- 
ized SQUID geometries that are amenable to analytical 
treatment. First, we consider a thin wire of diameter D 
in a circular loop of radius R with i? 3> D. We find that 
the mean-square flux induced in the SQUID by the TLS 
defects is given by 



2p\ 2 R 



(4) 



where a is the density of TLS surface defects on the su- 
perconducting wire. A factor 1/3 arises from a random 
angular distribution of the TLS magnetic moments. For 
the more realistic geometry of a thin-film superconduc- 
tor of width W and thickness 6 in a circular loop of ra- 
dius R with R ^S> W ^S> b, the surface currents J(x) at 
position —W/2 + A < x < W/2 — A are proportional to 
[1 — (2x/W) 2 ]^ 1 ^ 2 , where A is the penetration depth [15]. 
The currents fall away exponentially to zero at the edges 
±W/2. With the surface magnetic field being propor- 
tional to the surface current density, the mean-square 
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flux coupled to the SQUID is calculated as follows: 



<$ 2 ) = (w/6)p 2 p? B aR 



Jdx J 2 (x) 



2f4 



2 R 
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[Jdx J(x)Y 
ln(26H7A 2 ) 
2tt 



0.27 



(5) 
(6) 



The logarithmic term changes the prediction of Eq. Q 
by a factor ~ 1.8, with a reasonable fraction of the noise 
arising from fluctuators within a few penetration depths 
near the edges of the film. 

The major geometric dependence of the noise comes 
from the ratio R/W, the loop radius to width, with only 
a logarithmic dependence on the overall scale [16j. This 
feature of the model is compatible with the observation 
that the flux noise of /im-sized flux qubits is similar to 
that found for our 200 pm scale qubit, as the geometric 
ratio R/W is similar for these devices. 

The critical parameter determining the magnitude of 
the noise is the surface density of defect states a = pt, 
where t is the thickness of the surface oxide on the super- 
conducting film. The TLS defect density in amorphous 
oxide films can be extracted from measurements of the 
loss tangent of large-area tunnel junctions [I]. It is found 
that this defect density is compatible with bulk values 
obtained for a wide variety of amorphous oxides. We 
therefore take as an estimate of the TLS surface den- 
sity [4] pt = 1.0/^m 2 (/iGHz), twice that measured in 
tunnel junctions, to account for the thicker surface oxide 
t ~ 2 nm. To calculate the spectral density of the flux 
noise, we substitute S m t for p 2 B a in Eq. Using the pa- 
rameters R/W = 10 and T = 100 mK, we compute a flux 
noise spectral density S$(l Hz) = 1.1 x 10~ 3 (^<I>o) 2 /Hz, 
about 4 orders of magnitude smaller than the measured 
flux noise. 

Although we have not explicitly considered the noise 
contribution from other surfaces away from the super- 
conductor or dielectrics in crossover wiring, these small 
volumes cannot compensate for the large discrepancy be- 
tween the measured and calculated noise. The substrates 
are not likely candidates since they are typically crys- 
talline, and therefore have very low defect densities. In 
addition, defects at a Si/SiO x interface cannot account 
for the measured noise since our devices were made on 
sapphire substrates. 

Moreover, we note that the assumption that TLS fluc- 
tuation randomizes the defect magnetic moment is highly 
questionable because TLS defects in typical oxides are 
not considered to be magnetic. The above density of 
magnetic defect states is probably a gross overestimate, 
further exacerbating the discrepancy between the mea- 
sured flux noise and the noise calculated from TLS defect 
states. 

If spin noise is responsible for flux noise, we conclude 
that it must arise from a surface defect mechanism that 
is very different than that described by the standard TLS 



model, as it must have a defect /atomic-bond ratio that is 
about 10 4 times larger than for bulk TLS defects. Such 
a model would predict that specific heat measurements 
for amorphous materials would be dominated by surface 
states once the thickness of structures is less than about 
1 — 10 pm. Koch et al. have suggested [3] surface elec- 
tronic states as a possible candidate; unfortunately, the 
density of these defects has been estimated only at room 
temperature. A possible new mechanism has been pro- 
posed based on tunneling of conduction electrons into 
surface states [T7] . 

In conclusion, we have demonstrated a new measure- 
ment of 1// flux noise in superconducting qubits, which 
allows us to distinguish between flux and critical-current 
fluctuations. The magnitude of the measured noise is in 
good agreement with previous experiments, even though 
device parameters greatly differ. We have also theoreti- 
cally considered a spin-noise mechanism arising from fluc- 
tuating TLS. With the predicted magnitude in disagree- 
ment by over 4 orders of magnitude, we conclude that any 
model for spin noise must arise from a new mechanism 
based on a high density of defects. 
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